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Abstract 

Circinus X-1 exhibited a bright X-ray flare in late 2013. Follow-up observations with Chandra and 
XMM-Newton from 40 to 80 days after the flare reveal a bright X-ray light echo in the form of four well- 
defined rings with radii from 5 to 13 arcminutes, growing in radius with time. The large fluence of the 
flare and the large column density of Interstellar dust towards Circinus X-1 make this the largest and 
brightest set of rings from an X-ray light echo observed to date. By deconvolving the radial Intensity 
profile of the echo with the MAXI X-ray lightcurve of the flare we reconstruct the dust distribution 
towards Circinus X-1 into four distinct dust concentrations. By comparing the peak in scattering 
intensity with the peak intensity in CO maps of molecular clouds from the Mopra Southern Galactic 
Plane CO Survey we Identify the two innermost rings with clouds at radial velocity ^ —74kms“^ 
and ~ —81 kms“^, respectively. We Identify a prominent band of foreground photoelectric absorption 
with a lane of CO gas at ~ —32kms“^. From the association of the rings with individual CO clouds 
we determine the kinematic distance to Circinus X-1 to be Dcrx-i = 9-4+5 o kpc. This distance rules 
out earlier claims of a distance around 4kpc, implies that Circinus X-1 is a frequent super-Eddington 
source, and places a lower limit of F ^ 22 on the Lorentz factor and an upper limit of 0jet ^ 3° on 
the jet viewing angle. 

Subject headings: ISM: dust, extinction — stars: distances — stars: neutron — stars: individual 
(Circinus X-1) — techniques: radial velocities — X-rays: binaries 


1. INTRODUCTION 

I.l. X-ray Dust Scattering Light Echoes 

X-rays from bright point sources are affected by scat¬ 
tering off interstellar dust grains as they travel towards 
the observer. For sufficiently large dust column densi¬ 
ties, a significant fraction of the flux can be scattered 
into an arcminute-sized, soft X-ray halo. Dust scattering 
halos have long been used to study the properties of the 


properties of the X-ray source (e.g. 

Mathis & Lee 

1991 

Predehl & Schmitt||1995 

Corrales & Paereis||2U13|). 


delay experienced by the scattered X-ray photons rela¬ 
tive to the un-scattered X-rays can provide an even more 
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powerful probe of the ISM than a constant X-ray halo. 
In this case, the dust scattering halo will show temporal 
variations in intensity that can be cross-correlated with 
the light curve to constrain, for example, th e distance to 


the d ust clouds along the line-of-sight (e.g., Xiang et al. 


20111 . 

A^pecial situation arises when the X-ray source ex¬ 
hibits a strong and temporally well-defined flare (as ob¬ 
served, for example, in gamma-ray bursts and magne- 
tars). In this case, the scattered X-ray signal is not a 
temporal and spatial average of the source flux history 
and the Galactic dust distribution, but a well-defined 
light echo in the form of distinct rings of X-rays that 
increase in radius as time passes, given the longer and 
longer time delays (|Vianello et al.||2007|). 

The scattering geometry of X-rays reflected by an in¬ 
terstellar dust cloud is illustrated in Fig. For a source 
distance D and a layer of interstellar dust at distance 
xD from the observer (such that the fractional distance 
to the dust layer is a;), the scattered X-rays must travel 
an additional distance AD which depends on the ob¬ 


served off-a xis angle 6 (e.g. Mathis & Lee 1991 Xiang 

eTMlIMIj ): 


AD = xD (- l) + {1 - x) D (- l) 
\cos[0) ) \cos(a) ) 

^xDtfl X (\ — x)Do? xDO"^ 

2 ^ 2(l-a;) 

where we have used the small-angle approximation. The 
scattered signal will arrive with a time delay of 

, AD xDO^ 

At = -= 

c 


2c(l — x) 


( 2 ) 
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used throughout this paper. X-rays from a source at distance D 
scatter off a single dust layer at distance xD (such that the frac¬ 
tional distance is x). The observed angle of the scattered X-rays 
is while the true scattering angle is 0sc = 0 a. The scattered 
X-rays travel an additional path length of AD = ADi H- AD 2 . 



For a well-defined (short) flare of X-rays, the scattered 
flare signal propagates outward from the source in an 
annulus of angle 


2cAt(l — x) 


xD 


( 3 ) 


For a point source X-ray flux Fit) scattered by a thin 
sheet of dust of column density -/Vh, the X-ray intensity 
Iv of the dust scatter ing annulus at angle 9 observed at 
time tobs is given by (Mathis fc Lee|[T9M]) 


Iv{9, fobs) = -Xh 


-^I^(fobs At^) d(Jy 


[i-xY 


dfl 


( 4 ) 


where Atg is given by eq. (U and day/d^l is the differ¬ 
ential dust scattering cross section per Hydro gen atom, 
integrated over the gra in size distribution (e.g. Mathis & 
Lee|i991 Draine|2003 1. In the following, we will refer to 
the quantity JM}iday/dfl as the dust scattering depth. 

For a flare of finite duration scattering off a dust layer 
of finite thickness, the scattered signal will consist of a 
set of rings that reflect the different delay times for the 
different parts of the lightcurve relative to the time of the 
observation and the distribution of dust along the line of 
sight 

In tHe presence of multiple well-defined scattering lay¬ 
ers at different fractional distances x (with thickness 
^dust ^ D), each layer will generate a separate set of 
annuli. At a given time delay, the annuli of a fixed, short 
segment of the lightcurve are defined by the intersections 
of an ellipsoid of constant path length and each scatter¬ 
ing layer. The superposition of the annuli from all dust 
layers (clouds) and the different parts of the flare light 
curve will create a set of (partially overlapping) rings. 
This situation is illustrated in Fig. 

If the lightcurve of the X-ray flare that created the 
echo is known, an X-ray image of such a set of light- 
echo rings with sufficiently high angular resolution and 
signal-to-noise may be de-composed into contributions 
from a set of distinct dust layers at different distances. 
Such observations require both high angular resolution 
and sufficiently large collecting area and low background 


Because we will discuss delay times of several months, the 
much shorter duration of the X-ray observations may be neglected 
in the following discussion. 


Fig. 2. — Cartoon of the X-ray dust-scattering echo analyzed in 
this paper: A bright flare from a point source (center of the source 
plane) propagates towards the observer. Dust in dense interstellar 
clouds between observer and source scatters X-rays towards the 
observer, producing a light echo in the form of well defined rings. 
In this example, four clouds [a]-[d] shown in blue, green, yellow, 
and red, are located at different distances, illuminated by a bright 
flare (blue lightcurve) and producing four well defined rings in the 
image/source plane. 


to detect the rings. The Chandra X-ray observatory is 
ideally suited for such observations. 


1.2. Dust Echoes of Galactic X-ray Sources 

X-ray rings generated by dust scattering light echoes 
are extremely rare. To date, rings have been repo rted 
from five gamma-ray bursts (Vianello et al. 2007 and 
references therein). For Galactic X-ray sources, observa- 
tions of light echoes have been even more elusive. The 
soft gamma-ray repeater IE 1547.0-5408 was the first 
Galactic X-ray s ource to show a clea r set of bright X-ray 
light echo rings ( Tiengo et al.||2010 ). Two other sources 
have also been reported to show ligh t echo rings: the soft 
gamma-ray repeater SGR 1806-20 dSvirski e t al.||2011|) 
and t he fast X-ray transient IGR J17544-2619 (jiVlao et al. 
20141. However, the fluence of the initial X-ray hare 
in the latter two cases was too low to produce clearly 
resolved, arcminute-sca le rings and both detectio ns are 
tentative. In addition, McCollough et al.| (20131 found 
repeated scattering echoes from a single ISiVI cloud to¬ 
wards Gygnus X-3. 

In the case of GRBs, the source is known to be at cos¬ 
mological distances, which makes the scattering geome¬ 
try particularly simple, since the source is effectively at 
infinite distance and the scattering angle is equal to the 
observed ring angle on the sky, which simplihes eq. ^ 
to 


^GRBecho = 


2cAt 

hddust 


( 5 ) 


An observed ring radius ^GRBecho and delay time At 
therefore unambiguously determines the distance to the 
dust. 

However, in the case of Galactic sources, eq. .© for 
9{At, D,x) cannot be inverted to solve for the distance 
to either the dust or the source individually, unless the 
distance to either the source or the dust is known apriori. 
If a particular ring can be identified with a dust cloud 
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Fig. 3.— MAXI 2-4 keV lightcurve of Circinus X-1 at the time 
of the 2013 flare. Also shown are the median-times of our two 
Chandra and three XMM observations {XMM ObsID 0729560701 
was not used in this paper). Tq was chosen to correspond to the 
time of our first Chandra observation at (Tq =MJDi58oi= 56683). 
Red diamond s ymbols indicate periastron at orbital phase zero 
( |Nicolson|2007[ |. 

of known distance xD, the distance to the source can be 
determined from the known delay time At. 



15:22:00 21:00 20:00 

RA (J2000) 


1.3. Circinus X-1 

Circinus X-1 is a highly variable X-ray binary that has 
often been characterized as erratic. It has been difficult 
to classify in terms of canonical X-ray binary schemes, 


1 Stewart et al. 11991 

lOosterbroek et al. 1995 

Jonker et al. 

21)07 Calvelo et al. 

I2U12I). I'he presence of tv 

pe 1 X-ray 

bursts (I'i'ennant et al.|l986 

Linares et al.|2010 

) identifies 


a low magnetic held, supported by the presence of jets 
and the lack of X-ra y or radio pulses dStewart et 1993 


|2004a| 


Tudose et al.|2d06f Heinz et al. 2UU7 

Soleri et al.||2009| |l3eil et ai.||2010j). 

I'he transient nature and the low magnetic held sug¬ 
gested a classihcation of Circinus X-1 as an old, lo w-mass 
X-ray binar y with an evolved companion star (e.g. Shirey 
et al.|[T9^ ). This interpretation was in conhict with the 
observ ed orbital evolution time of order P/P ^ 3, 000 
years (Parkinson et al. 2003 Clarkson et al. 20041 to 


Fender et al. 


P/P ~ 20,000 years (Nicolson 2007) and the possible 
identihcation of an A-B type super-giant companion star 


(Jonker et al. 


2007), suggesting a much younger system 


age and a likely massive companion star. The orbital pe¬ 
riod of the s ystem is 16.5 days, with a likely eccentricity 
around 0.45 ( Jonker et al.||2007 ) 

The conhicting identihcations were resolved with the 
discovery of an arcminute scale su pernova remnant in 
both X-ray and radio observations (Heinz et al. 2013), 
placing an upper limit of r < 4, 600 years (Ucir/S^ kpc) 
on the age of the system, where ZJcir is the distance 
to the source. This upper limit makes Circinus X-1 the 
youngest known X-ray binary and an important test case 
for the study of both neutron star formation and orbital 
evolution in X-ray binaries. 

The distance to Circinus X-1 is highly uncertain. As a 
southern source deeply embedded in the Galactic disk, 
most standard methods of distance determination are 


Fig. 4.— Exposure-corrected image of Chandra ObsID 15801, 
smoothed with a FWHM=10.5” Gaussian. Red: 1-2 keV, green: 2- 
3 keV, blue: 3-5 keV. Point sources and the read streak produced by 
Circinus X-1 were removed. Four separate rings of dust scattering 
echoes are visible in this image. 


not available. It is also out of the range of VLBI par¬ 
allax measurements. Distance estimates to Circinus X-1 


range from 4kpc (laria et al. 2005) to llkpc (Stewart 


et al. 1991). A likely range ot L > = 8 — 10.5 tpc was 


proposed by Jonker & Nelemans (2004) based on the 
radius-expansion burst method, using properties of the 
observed type I X-ray bursts. Because important binary 
properties depend on the distance to the source (such 
as the Eddington fraction of the source in outburst), a 
more accurate determination of the distance is critical for 
a more complete understanding of this important source. 

Circinus X-1 is a highly variable source, spanning over 
four orders of magnitude in luminosity from its quiescent 
flux to its brightest flares. While the source was consis¬ 
tently bright during the 1990’s, it underwent a secular 
decline in flux starting in about 2000 and has spent the 
past 10 years below the detection thresholds of the Rossi 
X-ray Timing Explorer All Sky Monitor and the MAXI 
all sky monitor aboar d the International Space Station 
(Matsuoka et al. |2009l|^^ with the exception of sporadic 
flares, S( 


near-Eddington flares, sometimes in excess of one Crab 
in flux (throughout the paper we will refer to the Edding¬ 
ton limit of TEdd.NS ~ 1-8 x 10^® ergs s“^ for a 1.4 Mq 
neutron star). These flares typically occur within indi¬ 
vidual binary orbits and are characterized by very rapid 
rises in flux near periastron of the binary orbit and often 
rapid flux decays at the end of the orbit. For example, 
during the five years it has been monitored by MAXI, 
Circinus X-1 exhibited ten binary orbits with peak 2- 


The ASM was one of three instruments aboard the Rossi 
X-ray Timing Explorer, MAXI is an experiment aboard the Inter¬ 
national Space Station-, both have only one mode of operation. 
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15:22:00 21:00 20:00 19:00 

RA (J2000) 


Fig. 5.— Exposure-corrected and background-subtracted image 
of XMM-Newton ObsID 0729560501 (combining EPIC MOS and 
PN data), smoothed with a FWHM=35” Gaussian. Colors rep¬ 
resent the same energy bands as in Fig. Point sources have 
been removed. The over-plotted white-dasned line corresponds to 
the 15K kms“^ contour of the —33.6 to —29.6 kms“^ image 

shown in Fig. |19| 


10 keV fluxes in excess of ^ 10“® ergs cm“^ s 
orbits with peak flux at or above one Crab. 

The large dynamic range of the source flux, the well- 
defined duration of the X-ray flares, and the low mean 
flux during the quiescent periods over the past decades 
make Circinus X-1 an ideal source for observations of X- 
ray light echoes. The source is located in the Galactic 
plane, at Galactic coordinates I = 322.12°, 6 = 0.04°. 
The large neutral hydrogen column density of Vh ~ 2 x 
10^^ cm“^, inferred from photo-electric absorp tion of the 
X-ray spectr a of both the neutron star (e. g. [Brandt fc 
Schul z 20001 and the supernova remnant (Heinz et al. 
2013|, further increases the likelihood to observe light 
echoes. 

In late 2013, Circinus X-1 underwent an extremely 
bright flare, preceded and followed by long periods of 
very low X-ray flux. In this paper, we report the detec¬ 
tion of the brightest and largest set of X-ray dust scat¬ 
tering rings observed to date, generated by the light echo 
from the 2013 flare. The spatial variation of the rings al¬ 
lows us to identify the interstellar clouds responsible for 
the different rings, which determines the distance to the 
X-ray source, Circinus X-1, to roughly 10% accuracy. 

We will describe the observations and data analysis in 
^ Section describes the analysis of the X-ray, CO, 
and HI data and the procedure for deconvolving the X- 
ray signal into the line-of-sight dust distribution. In 0 
we derive a new distance to Circinus X-1 and discuss the 
consequences for the physics of this unique X-ray binary 
and for the physics of dust scattering. Section briefly 
summarizes our findings. 


2. OBSERVATIONS AND DATA REDUCTION 



Fig. 6. — Un-absorbed 2-5 keV scattering intensity in units of 
10“^ ergscm“^ s“^ Sr“^ across the C/iandra field-of-view, derived 
from spectral fits to Chandra ObsID 15801 on a rectangular grid of 
8x8 squares per CCD, labeled by row and column from 1 through 
24 and A through P, respectively. Photo-electric absorption was 
removed to isolate the dust-scattered intensity. 


2.1. X-Ray Observations 

Circinus X-1 was observed with Chandra, XMM- 
Newton, and Swift in early 2014 after a bright X-ray 
flare of the source in late 2013. Figure shows the time 
of the Chandra and XMM-Newton observations overlaid 
on the MAXI lightcurve of the source. The source ex¬ 
hibited a bright flare during three binary orbits approx¬ 
imately 50_to_JT0du;TO before the first Chandra obser¬ 
vation (Asai et al. 20141, with a quiescent binary orbit 
interspersed during the first half of the flare. The source 
flare reached fluxes in excess of one Crab, with a total 2-4 
keV fluence of ^ 2-4 = 0.025 ergs cm“^. The lightcurve of 
the flare has three characteristic peaks and a pronounced 
gap. This is the most energetic flare Circinus X-1 exhib¬ 
ited since it entered its long-term quiescent state in 2006, 
and is comparable in peak flux to the brightest flares the 
source has undergone even during its long-term outburst. 

Before and during the Chandra campaign. Swift ob¬ 
served the source twelve times for a total exposure of 
12 ksec in order to safely schedule the Chandra observa¬ 
tions. The history of bright flares by Circinus X-1, which 
can exceed the ACIS dose limit even during moderate ex¬ 
posures presents a significant scheduling challenge, and 
frequent monitoring of the source with Swift was used 
to confirm that the source was safely below the max¬ 
imum flux allowable under Chandra instrument safety 
rules. The short total exposure time and the fact that 
the exposures were spread over a time window of 11 days 
make the Swift data unsuitable for the analysis presented 
in this paper. 


2.1.1. Chandra Observations 

In the 3.5 binary orbits preceding the first Chandra 
observation, as well as during the months prior to the 
flare, the source was quiescent, below the MAXI detec- 
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tion threshold. We observed the source on two occasions 
with Chandra, on January 25, 2014 for 125 ksec, and on 
January 31, 2014 for 55 ksec, listed as ObsID 15801 and 
16578, respectively. The point source was at very low 
flux levels during both observations. Circinus X-1 was 
placed on the ACTS S3 chip. 

Data were pipeline processed using CIAD software ver¬ 
sion 4.6.2. Point sources wer e identified using the 
wavdetect (Freeman et al.||2002) task and ObsID 16578 
was reprojected to match the astrometry of ObsID 15801. 
For comparison and analysis purposes, we also repro¬ 
cessed ObsID 10062 (2009) with CIAO 4.6.2 and repro¬ 
jected it to match the astrometry of ObsID 15801. 

We prepared blank backgroun d images following the 


standard CIAD thread (see also Hickox & Markevitch 


20061, matching the hard (>10 keV) X-ray spectrum of 
the background hie for each chip. 

Figure shows an exposure-corrected, background- 
subtractea three-color image of the full ACIS held-of 
view captured by the observation (ACIS chips 2,3,6,7, 
and 8 were active during the observation), where red, 
green, and blue correspond to the 1-2, 2-3, and 3-5 keV 
bands, respectively. Point sources were identihed using 
wavdetect in ObsIDs 10062, 15801, and 16578, source 
lists were merged and point sources were removed from 
imaging and spectral analysis. The read streak was re¬ 
moved and the image was smoothed with a 10.5” full- 
width-half-max (FWHM) Gaussian in all three bands. 

The image shows the X-ray binary point source, the X- 
ray jets (both over-exposed in the center of the image), 
and the supernova remnant in the central part of the im¬ 
age around the source position at 15:20:40.9,-57:10:00.1 
(J2000). 

The image also clearly shows at least three bright rings 
that are concentric on the point source. The first ring 
spans from 4.2 to 5.7 arcminutes in radius, the second 
ring from 6.1 to 8.2 arcminutes, and the third from 8.3 
to 11.4 arcminutes in radius, predominantly covered by 
the Eastern chips 2 and 3. We will refer to these rings 
as rings [a], [b], and [c] from the inside out, respectively. 
An additional ring-like excess is visible at approximately 
13 arcminutes in radius, which we will refer to as ring 

As we will lay out in detail in ® we interpret these 
rings as the dust-scattering echo from the bright flare 
in Oct.-Dec. of 2013, with each ring corresponding to 
a distinct concentration of dust along the line-of-sight 
to Circinus X-1. A continuous dust distribution would 
not produce the distinct, sharp set of rings observed by 
Chandra. 

The rings are also clearly visible in ObsID 16578, de¬ 
spite the shorter exposure and the resulting lower signal- 
to-noise. In ObsID 16578, the rings appear at ~4% larger 
radii, consistent with the expecta tion of a dust ech o mov¬ 
ing outward in radius [see also [3.2 and eq. ®]- The 
rings are easily discernible by eye m the energy bands 
from 1 to 5 keV. 

Even though the outer rings are not fully covered by 
the Chandra field of view (abbreviated as FOV in the 
following), it is clear that the rings are not uniform in 
brightness as a function of azimuthal angle. There are 
clear intensity peaks at [15:21:00,-57:06:30] in the inner 
ring [a] and at [15:20:20,-57:16:00] in ring [b] of ObsID 
15801. Generally, ring [c] appears brighter on theSouth- 



Fig. 7.— Distribution of photo-electric absorption column A^h 
derived from the spectral fits shown in Fig. [^ A clear excess of 
absorption is found in the lower half of the image. Over-plotted is 
thel5Kkms“^ contour of the —33.6 to—29.6 kms“^ image 

shown in Fig. |19| 


Eastern side of the image. 

The deviation from axi-symmetry observed here dif¬ 
fers from almost all other observations of dust scattering 
signatures, which typically appear to be very uniform 
in azimuth. For example, a detailed investigation of the 
dust scattering halo of Cygnus X-2 found that th e profile 


k, Smith 

20131 

nature o 

non-^ 


20131. The only comparable observational sig- 


echo off a bok glo bule observed for Cygnus X-3 (jMcCol- 
lough et ^[2013 1 


In addition to the light echo, a blue ’’lane” is visible 
across the lower part of the image. As we will show be¬ 
low, this feature is the result of photoelectric foreground 
absorption by an additional layer of dust and gas closer 
to the observer than the clouds producing the light echo. 

2.1.2. XMM-Newton Observations 

XMM-Newton observed the source on 17 February, 
2014 (ObsID 0729560501), 4 March, 2014 (ObsID 
0729560601), and 16 March, 2014 (ObsID 0729560701) 
for 30 ksec each. ObsID 0729560501 was not affected 
by any background flares. The second and third ob¬ 
servation were strongly affected by background particle 
flares, eliminating the use of PN data completely for 
ObsID 0729560601, and the effective exposure for the 
MOS detectors after flare rejection was cut to approxi¬ 
mately 15ksec in both cases. The resulting data are noisy, 
suggesting that additional background contamination is 
present, making the data not usable for detailed, quan¬ 
titative image processing. We restrict the discussion in 
this paper to data set 0729560501 (with the exception of 
Fig.pf. 

MOs and PN data were reduced using SAS software 
version 13.5. Data were prepared using the ESAS package, 
following the standard diffuse source analysis guidelines 
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Fig. 8.— Un-absorbed excess 2-5 keV scattering intensity in units 
of ergscm“^ s“^ Sr“^ across the XMM field-of-view, derived from 
spectral fits in annular sections identical to those plotted in Fig.f^ 
Photo-electric absorption was removed to isolate the dust-scattered 
intensity. Clear peaks in the surface brightness are located ring 
sections section [all] and [b6]. 


from the Diffuse Analysis CookboolPI addition to 
point sources identified by the ESAS cheese task, point 
sources identified in the three Chandra exposures (10062, 
15801, 16578) were removed before image and spectral 
extraction. The energy range from 1.4-1.6 keV around 
the prominent instrumental Aluminum line at 1.5 keV 
was removed from image processing, while spectral fits 
determined that the contribution from the Si line at 1.75 
keV was sufficiently small to ignore in image processing. 
Proton and solar-wind charge exchange maps were cre¬ 
ated using spectral fits to the outer regions of the image 
and subtracted during image processing. 

Figure [^shows a three color exposure-corrected im¬ 
age of OMD 0729560501 in the same bands as Fig. 
smoothed with a 35” FWHM Gaussian. Despite the 
lower spatial resolution of XMM-Newton, the image 
clearly shows the three rings identified in the Chandra 
images. The median radii of the rings have increased 
by ^16% relative to ObsID 15801, corresponding to the 
longer time delay between the flare and the XMM obser¬ 
vation. 

Overlaid on the figure is a grid of three annuli covering 
the rings, broken into twelve angular sections of 30° each. 
We will use this grid in the discussion below to identify 
different areas in the image and in corresponding CO 
maps of the FOV. The rings are labeled [a],[b], and [c] 
from inside out, and t he se ctions are labeled 1-12 in a 
clock-wise fashion. See §3.1| and Fig.for further details 
on the grid and the nomenclature used in this paper. 

2.2. CO Data 

We used CO data from the Mopra Southern Galac¬ 
tic Plane CO Survey, described in Burton et al. 2013, 
with the particular data set used here, for the Circi- 
nus X-1 region (1=322), coming from the forthcoming 
public data release in Braiding et al., 2015 (in prepara- 

ftp://xmm.esac.esa.int/pub/xmm-esas/xmm-esas.pdf 



Fig. 9.— Distribution of photo-electric absorption column TVh 
derived from the spectral fits shown in Fig. A clear excess of 
absorption is found in the lower half of the image, with a peak at 
position [c5]. 


tion). In this data set emission from the J=l-0 line of the 
three principal isotopologues of the CO molecule (^^CO, 
i^CO, Ci®0, at 115.27, 110.20 and 109.78 GHz, respec¬ 
tively) has been measured at 0.6 arcmin, 0.1kms“^ res¬ 
olution with the 22m diameter Mopra millimetre wave 
telescope, located by Siding Spring Observatory in New 
South Wales in Australia. 

By continuously scanning the telescope while making 
the observations (on-the-fly mapping) large areas of sky 
can be readily mapped by Mopra. Data cubes covering 
1° X 1° of sky typically take 4 nights with the telescope to 
obtain, and achieve a 1 sigma sensitivity of about 1.3K 
and 0.5K per 0.1kms“^ velocity channel for the 12CO 
and the (13CO, C180) lines, respectively. Such a data 
cube, covering the G322 region {I = 322° — 323°, b = 
—0.5° to -1-0.5°), was used for the analysis presented here. 


2.3. 21 cm Data 

We use publicly available 21 cm radio data from 
the Parkes/ATCA Southern Galactic Plane survey 
(McClure-Griffiths et al. 2005) to estimate the column 
density of neutral atomic hydrogen towards Circinus X-1. 
We use CO data to identify velocity components because 
the velocity dispersion of the HI velocity components is 
too large to unambiguously resolve individual clouds. 


3. ANALYSIS 
3.1. X-ray Spectral Fits 

We performed spectral fits of data from Chandra Ob¬ 
sID 15801 and XMM ObsID 0729560501, with results 
from the spectral analysis presented in Figs. Be¬ 

cause the XMM-Newton FOV is contiguous andcovers 
the rings entirely, we extracted spectra for each ring of 
XMM ObsID 0729560501 in annuli spanning the radial 
ranges [4.7’-7.1’,7.F-9.7’,9.7’-13.7’], with each ring di¬ 
vided into twelve sections. Rings are labeled [a] through 
[c] (with lower case letters denoting XMM rings), and 
the layout of the grid is shown in Fig. We divided the 
Chandra FOV into a rectangular grid7 with each CCD 
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covered by an 8x8 grid of square apertures. The layout of 
the grid on the ACIS focal plane is shown in Fig. with 
columns running from [A] through [P] and rows from 1 
to 24 (with capital letters denoting Chandra columns). 

Spectra for Chandra ObsID 15801 were generated us¬ 
ing the specextract script after point source removal. 
Background spectra were generated from the blank sky 
background files. No additional background model was 
required in fitting the Chandra data, given the accurate 
representation of the background in the black sky data. 

Spectra for XMM-Newton ObsID 0729560501 were ex¬ 
tracted from the AMM-MOS data using the ESAS pack¬ 
age and following the Diffuse Analysis Cookbook. In 
spectral fitting of the XMM data, we modeled instrumen¬ 
tal features (A1 and Si) and solar-wind charge exchange 
emission by Gaussians with line energies and intensities 
tied across all spectra and line widths frozen below the 
spectral resolution limit of AMM-MOS at cr = 1 eV. We 
modeled the sky background emission as a sum of an 
absorbed APEC model and an absorbed powerlaw. Inten¬ 
sities of both background components were tied across 
all ring sections. 

We modeled the dust scattering emission in each re¬ 
gion as an absorbed powerlaw, with floating absorption 
column and normalization, but with powerlaw index tied 
across all regions, with the exception of the tiles cover¬ 
ing the central source and the supernova remnant in the 
Chandra image, which span a square aperture from posi- 
ton [J14] to [M18]. Tying the powerlaw photon index P 
across all spectra is appropriate because (a) the scatter¬ 
ing angle varies only moderately across the rings, with 


Oar. = 


2cAt 


(1 — x) xD 


( 6 ) 


and (b) for the scattering angles considered here {0sc > 
10'), the energy dependence of the scattering cross sec- 
tion is indep endent of dsc, with da/dCl oc E~'^ (e.g. 
Draine |2003|). While the source spectrum might have 
varied during the outburst, we are integrating over suffi¬ 
ciently large apertures to average over possible spectral 
variations. 

While the source spectrum during outburst (the echo 
of which we are observing in the rings) was likely more 
complicated than a simple powerlaw, we have no direct 
measurement of the detailed flare spectrum and the qual¬ 
ity of the spectra does not warrant additional complica¬ 
tions, such as the addition of emission lines which are 
often present in the sp ectra when the source is bright 
( Brandt fc Schulz||2000 ). Dust scattering is expected to 
steepen the source spectrum by AF = 2, and the mea¬ 
sured powerlaw index of 


r = 4.00 ±0.03 


( 7 ) 


for the scattered emission is consistent with the typically 
soft spectrum of F ~ 2 the source displays when in out¬ 
burst. 

The spectral fits resulting from the model are sta¬ 
tistically satisfactory, with a reduced chisquare of 
X?ed, 0729560501 = 839.68/807 d.o.f. = 1.04 for the com¬ 
bined fit of the entire set of spectra for ObsID 0729560501 
and Xred 15801 = 839.68/807 d.o.f. = 1.05 for ObsID 



Fig. 10.— Radial X-ray intensity profile of Chandra ObsID 15801 
in the 2-5keV band (chosen to roughly match the MAXI 2-4 keV 
band), shown in black. The intensity profile of Chandra ObsID 
10062 was subtracted to remove the supernova remnant and resid¬ 
ual quiescent background. Overlaid is the full reconstructed radial 
intensity profile from the deconvolution (gray), as well as the con¬ 
tributions from the four dist inct dust concentrations along the line 
of sight identified in Fig. |12| and the powerlaw model of the diffuse 
instantaneous dust scattering halo (purple line). The hatched area 
shows the region where the deconvolution becomes unreliable be¬ 
cause part of the kernel falls outside of the Chandra FOV or onto 
the supernova remnant. Top right insert: Dust scattering kernel 
A"(^/0med) [see eq. J^] used for the deconvolution, derived from 
the MAXI 2-4 keV ^htcurve during the flare shown in Fig. 
normalized by the total fiuence of the outburst. 

I580I. We display the unabsorbe^^ 2-5keV inten¬ 
sity maps of Chandra ObsID I580I and XMM ObsID 
0729560501 in Figs. and respectively. The unab¬ 
sorbed intensity provides a direct measure of the scatter¬ 
ing intensity (unaffected by effects of foreground absorp¬ 
tion) and is directly proportional to the dust scattering 
depth. In Figs. and we show maps of the photo¬ 
electric absorption column density Nn determined from 
the spectral hts for Chandra ObsID 15801 and AMM Ob¬ 
sID 0729560501, respectively. 

Four properties of the spectral maps stand out that 
confirm the qualitative description of the color images in 
Figs. |4] and [5| 

1. A clear intensity peak at position [all] in Fig.j^and 
[119] in Fig. [^corresponds to the brightness peaks 
seen at that position in Figs. I^and The intensity 
is a factor of two larrer than in the neighboring 
section [al2] in Fig. [^ with a similar increment 
relative to the neighboring sections in Fig. 

2. There is a clear intensity peak at position [b5] 
and [b6], with [b6] being approximately two times 
brighter than neighboring [b7] in Fig. ^ Ring [b] 
is brighter on the Western side of the aMM FOV. 
This peak of ring [b] corresponds to the brightness 
peak in tiles [L9], [NIO], and [OlO] in Fig. 

3. The Southern section of ring [c] is brighter than 
the Northern section, with the brightest scattering 
emission in sections [c5:c8] in Fig.® Ring [d] over¬ 
laps significantly with ring [c], and^e color images 

Unabsorbed intensities are calculated from the model fits by 
removing the effects of photo-electric absorption before determin¬ 
ing model fluxes. 










































and intensity maps cannot distinguish the angular 
distribution of the two rings. 

4. There is a clear excess of foreground absorption in 
the Southern sections of the FOV, running from 
tile [P7] to [B 6 ] in Fig. Figure 1^ shows the 
strongest absorption in sections [c5], [d 6] and [b7]. 
This corresponds to the blue absorption lane that 
runs across the color images in Figs, [i] and 

The rings clearly deviate from axial symmetry about 
Circinus X-1, which requires a strong variation in the 
scattering dust column density across the image. Be¬ 
cause the rings are well-defined, sharp features, each ring 
requires the presence of a spatially well separated, thin 
distribution of dust (where thin means that the extent 
of the dust sheet responsible for each ring is significantly 
smaller than the distance to the source, by at least an 
order of magnitude). 

We present a quantitative analysis of the dust distri¬ 
bution derived from the rings in |3.2[ but from simple 
consideration of the images and spectral maps alone, it 
is clear that the different dust clouds responsible for the 
rings have different spatial distributions and cannot be 
uniform across the XMM and Chandra FOVs. In partic¬ 
ular, the peak of the cloud that produced the innermost 
ring (ring [a]) must have a strong column density peak 
in section [all], while the cloud producing ring [b] must 
show a strong peak in sections [b5] and [b 6 ]. The ex¬ 
cess of photo-electric absorption in the lower half of the 
image, with a lane of absorption running through the 
sections [c5], [b 6 ], and [b7], requires a significant amount 
of foreground dust and gas (an excess column of approx¬ 
imately A^h ~ 10 ^^ cm“^ relative to the rest of the field). 

3.2. X-ray Intensity Profiles and Dust Distributions 
3.2.1. Radial intensity profiles 

In order to determine quantitatively the dust distribu¬ 
tion along the line of sight, we extracted radial intensity 
profiles of Chandra ObsID 10062, 15801, and 16578 in 
600 logarithmic radial bins from 1 ’ to 18’.5 in the energy 
bands 1-2 keV, 2-3 keV, and 3-5 keV (after background 
subtraction and point-source/read-streak removal). We 
extracted the radial profile of XMM ObsID 0729560501 
in 300 logarithmic bins over the same range in angles. 
We subtracted the quiescent radial intensity profile of 
ObsID 10062 to remove the emission by the supernova 
remnant and residual background emission. 

Because of the variation in photo-electric absorption 
across the field, we restricted analysis to the 2-5 keV 
band, which is only moderately affected by absorption. 
We plot the resultin g ra dial intensity pro file of Chandra 
ObsID 15801 in Fig. The plot in Fig. [T0| shows four 
echo rings labeled [a] tmough [d] which correspond to the 
visually identified rings in the previous sections. Several 
bright sharp peaks in / 2 - 5 kev(^) are visible that can also 
be seen as distinct sharp rings in Fig. suggesting very 
localized dust concentrations. For comparison, the ra dial 
profile of XMM ObsID 0729560501 is plotted in Fig. [TT| 

3.2.2. Modeling and removing the diffuse instantaneous dust 
scattering halo 

In general, the dust scattering ring echo will be super¬ 
imposed on the lower-level instantaneous diffuse dust 



Fig. 11.— Radial X-ray intensity profile of XMM ObsID 
0729560501 in the 2-5 keV band, following the same labeling and 
convention as Fig. |10| The purple line shows the best fit powerlaw 
profile for the instantaneous dust scattering halo (removed before 
the deconvolution), the colored lines show the contributions from 
rings [a]-[d] derived from the radial deconvolution of the profile. 
Overplotted as dashed lines are the radial profiles from rings [a]-[d] 
predicted by the deconvolution of the Chandra ObsID 15801 radial 
profile. Top right insert: Dust scattering kernel R’(0/^med) [see 
eq. 01 used for the deconvolution. 


1.8 X 10 ^^ergss ^ cm ^ for ObsID 16578 
1.3 X 10“^^ ergss“^ cm“^ for 0729560501, with 


scattering halo produced by the (weak) emission of the 
point source during and just prior to the observation, 
which should be removed before analysis of the echo. The 
2-5 keV point source flux during and before ObsID 15801 
was very low (F2_5keV.i580i < 6 x 10 “^^ ergs“^ cm”^)^ 
resulting in essentially negligible halo emission, while the 
two later observations showed larger point source flux 
{F2-5keV 
and F ~ 

correspondingly brighter diffuse halo emission). 

We model the instantaneous dust scattering halo as a 
powerlaw in ring angle, /halo cx: 9^, appropriate for the 
large halo angles under consideration and consistent with 
the steady dust scattering halo prof ile for Circinus X - 1 
derived from Chandra ObsID 706 in Heinz et al. (20131. 

We determine the powerlaw index to be 77 ~ —1.5 by 
fitting the intensity profile of ObsID 0729560501 in the 
radial range between the edge of the supernova remnant 
at 9 > 2.5' and the inner edge of ring [a] at 9 < 5'. This 
value of t] is consistent w ith the reference du st scattering 
halo of ObsID 706 from Heinz et al. (2013), which has 
an index of 77 ~ —1.55 in the 2-5 keV band. 

We determine the powerlaw normalization for ObsID 
15801 and 16578 by fits to each radial profile between the 
outer edge of the remnant at 9 > 2.5' and the inner edge 
of ring [a] at 0 < 4'. We overplot the powerlaw fits to 
the instantaneous halo as purple lines in Fig. [T 0 ]and[TT 
We subtract the halo e missi on be fore further analysis of 
the ring echo in §3.2 4| and §3.2.5[ 


3.2.3. Deriving the dust distribution 

Given the well-sampled MAXI lightcurve of the 2013 
flare, we can construct the light echo intensity profile 
that would be produced by a thin sheet of dust at a 
given relative dust distance x and observing time fobs 
from eqs. (1^ and (|^. In general, the different parts of 
the flare will generate emission at different radii (echoes 
of earlier emission will appear at larger angles 9). 
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The scattering angle and thus the scattering cross sec¬ 
tion will vary for different parts of the lightcurve, which 
must be modeled. Given the large scattering angles of 
9sc ^ 1400" of the rings and photon energies above 2 
keV, it is safe to assume that the cross section is in the 
large angle lim it where da/dfl oc oc with a ~ 4 


(Draine 2003), which we will use in the following. The 
exact value of a depends on the unknown grain size dis¬ 
tribution. By repeating the analysis process with differ¬ 
ent values of a, ranging from 3 to 5, we have verified that 
our results (in particular, the location of the dust clouds 
and the distance to Circinus X-1) are not sensitive to the 
exact value of a. 

For a thin scattering sheet at fractional distance x, the 
median ring angle 0„ied (ke., the observed angle 9 for the 
median time delay Atmed of the flare) is 


^med — 




xD 


( 8 ) 


where Aimed is a constant set only by the date of the 
observation and the date range of the flare. The invariant 
intensity profile Al 2 - 5 keV produced by that sheet depends 
on the observed angle 9 only through the dimensionless 
ratio z = 9/9med' 


K. 


2 —5keV 


(z,x) 


= TMAXl(Ai„ied2:^) 


Ah,x da 

(l-a:)2 dn 


,-4 


(9) 


s=e„ 


We can therefore decompose the observed radial profile 
into a series of these contributions from scattering sheets 
along the line of sight, each of the form K2-5kev{S/9nied)- 
Because we chose logarithmic radial bins for the intensity 
profile, we can achieve this by a simple deconvolution of 
l 2 - 5 keV with the kernel K 2 - 5 keV- The inset in Fig. [TO" 
shows K{9/9^ed) for a = 4. 


3.2.4. Deconvolution of the Chandra radial profiles 

To derive the dust distribution that gave rise to the 
light echo, we performed a maximum likelihood decon¬ 
volution of the Chandra radial profiles with this kernel 
into a set of 600 dust scattering screens along the line 
of sight. We used the Lucy-Richardson deconvolution al¬ 
gorithm (Richardson 1972j_ Lucy 1974) implemented in 
the IDL ASTROLIB library*^] In the region of interest 
(where the scattering profile is not affected by the chip 
edge and the supernova remnant), the radial bins are 
larger than the 50% enclosed energy width of the Chan¬ 
dra point spread function (PSF), such that the emission 
inside each radial bin is resolvecpT| 

Figure shows the reconstructed intensity profile 
from this distribution for ObsID 15801 (derived by re¬ 
convolving the deconvolution with 7^2-5 keV and adding 
it to the powerlaw profile for the instantaneous dust scat¬ 
tering halo), which reproduces the observed intensity 
profile very well, with only minor deviations primarily 


http://idlastro.gsfc.nasa.gov/ftp/pro/image/max_likelihood.pro 
We verified that our results are robust against PSF effects by 
repeating the deconvolution after smoothing both the profile and 
the kernel with Gaussians of equal width; this mimicks the effects 
PSF smearing would have on the deconvolution. The smoothing 
does not affect our results significantly even for widths an order of 
magnitude larger than the Chandra PSF. 


Radial Velocity V^sr I km .s '| 
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Fig. 12.— Maximum likelihood d econ volution of the 2-5 keV ra¬ 
dial intensity profile shown in Fig. |10| plotted as a histogram of 
scattering depth along the line of sight (rebinned onto a uniform 
grid in 9). Shown in color (left-hand y-axis) is the invariant effec¬ 
tive scattering depth (1 — x)~^Nuda/dfl, where x is the fractional 
distance of the dust sheet relative to the distance to Circinus X-1 
and da/dti is the scattering cross section evaluated at the mean 
scattering angle of the sheet and at 2 keV. The distribution shows 
four main components (labeled [a]-[d] in colors blue, green, yel¬ 
low, and red, respectively). Also shown in gray is the inferred 
absolute scattering depth N^da/dQ, for the best fit distance of 
A*cirX—1 = 9.4kpc (right-hand y-axis), along with the absolute 
distance to the scatterer and the associated radial velocity due to 
Galactic rotation (top axes in gray). 


for rings [a] from 4’-6’, while the profile outward of 6’ is 
matched to w ithin the statistical uncertainties. 

In Fig. |12[ we plot the inferred dust distribution along 
the line ofaght as a function of median ring angle 0med, 
derived from our deconvolution of the intensity profile 
of ObsID 15801, and rebinned onto a uniform grid in 
9. The plotted quantity is a histogram of the invariant 
differential scattering depth AH(l—a;)^ da/dU,\g^g^^^ per 
bin, which is independent of the distance to Circinus X-1. 

The dust distribution clearly shows four main dust con¬ 
centrations toward Circinus X-1, which we plotted in dif¬ 
ferent colors for clarity. We fit each dust component with 
a Gaussian to determine the mean location and width, 
with fit values listed in Table Given the longer median 
time delay for ObsID 16578, the rings should be observed 
at larger median angles by a factor of 

^^med,16578 _ / Aimed,16578 _ /6.084 X 10® S _ ^ 

^^med,15801 y At^ed,15801 V 5.628 X 10® S 

( 10 ) 

or 4%, which is consistent with the observed increase in 
^ring.med, within the typical width of each cloud. 

We label the identified dust clouds as components [a], 
[b], [c], and [d] (plotted in green, blue, yellow, and red, re¬ 
spectively) according to the naming convention adopted 
above, and identify them with the main X-ray rings seen 
in the images. To show that this identification is appro¬ 
priate, we overplot the contribution of each dust cloud 
to the total intensity profile in the associated color in 
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TABLE 1 

Gaussian fits to dust clouds identified in the radial 

DECONVOLUTION IN FiG. fT^ 


Ring: 

^med 

O-Q 

do- ro^—lla 



Chandra ObsID 15801 

[a] 

5.37' ± 0.02' 

0.21' ± 0.02' 

52.56 ± 7.62 

0.83 

[b] 

7.63' ± 0.01' 

0.19' ± 0.01' 

47.86 ±4.77 

0.70 

[c] 

10.64' ± 0.01' 

' 0.22' ± 0.01' 

42.09 ±3.62 

0..' 

[d] 

12.85' ± 0.04' 

' 0.34' ± 0.04' 

19.60 ±3.17 

0.4." 

Chandra ObsID 16578 

[a] 

5.50' ± 0.02' 

0.15' ± 0.02' 

30.25 ±6.03 

0.83 

[b] 

7.91' ± 0.02' 

0.20' ± 0.02' 

36.16 ±4.18 

0.70 

[c] 

11.10' ± 0.02' 

' 0.26' ± 0.02' 

39.43 ±3.29 

0.55 

[d] 

13.49' ± 0.02' 

' 0.10' ± 0.02' 

6.51 ± 1.68 

0.4.". 

XMM ObsID 0729560501 

[a] 

6.00' ± 0.05' 

0.90' ± 0.06' 

29.31 ± 2.32 

0.84 

[b] 

8.72' ± 0.02' 

0.52' ± 0.02' 

25.56 ± 1.06 

0.71 

[c] 

12.04' ± 0.02' 

' 0.73' ± 0.02' 

23.54 ± 0.92 

0.56 


Note. — Cloud centroids correspond to the median time 
delay of Atisgoi = 5.628 X 10® s, Af igsyg = 6.203 x 10® s, 
and At072956050i = 7.925 x 10® s for ObsID 15801, 16578, and 
0729560501, respectively. We do not include ring [d] for XMM Ob¬ 
sID 0729560501 because its centroid falls outside of the MOS FOV. 
^ Invariant scattering depth, evaluated at the median scattering 
angle and at 2 keV energy. 

® relative distance x to the dust cloud responsi ble for the respective 
ring using the best fit distance of 9.4 kpc [eq. 


Fig. 10 Rings [a] and [b] are produced primarily by 
scattering from clouds [a] and [b], while ring [c] and [d] 
overlap significantly and contain contributions from both 
clouds. 

We tested the uniqueness of the deconvolution proce¬ 
dure using synthetic radial profiles and verified that the 
procedure faithfully reproduces the dust profile for ran¬ 
domly placed clouds down to a mean flux of about 50% 
of the background, for counts rates similar to the ones 
observed in Chandra ObsID 15801. Thus, while we can¬ 
not rule out the presence of additional weak dust rings 
below the detection threshold, the identification of the 
four main dust concentrations producing rings [a]-[d] is 
robust. 


3.2.5. Deconvolution of the XMM radial profile 

The XMM data are noisier due to (a) the shorter ex¬ 
posure, (b) the lower scattering intensity given the in¬ 
creased scattering angle at the correspondingly longer 
time delay, (c) the larger background intensity that 
is more difficult to model than the stable Chandra 
background, and (d) the brighter dust scattering halo 
due to the increased point source flux during ObsID 
0729560501. 

Despite the lower spatial resolution and significantly 
lower signal-to-noise ratio of the data, the rings are still 
discernible in the radial profile in Fig. We decon¬ 
volved the residual intensity profile following the same 
procedure used for the Chandra data While the 300 log¬ 
arithmic radial bins used for the deconvolution are nar¬ 
rower than the 50% enclosed energy width of XMM, we 
verified that the deconvolution was robust by confirming 
that a strongly smoothed radial intensity profile and/or a 


smoothed kernel lead to consistent locations and scatter¬ 
ing depths within the measurement uncertainties for the 
dust screens responsible for rings [a]-[c].. The resulting 
contributio ns o f rings [a] through [d] are overplotte d in 
color in Fig.f^ following the same convention as Fig.fTO} 
Gaussian fits to rings [a], [b], and [c] are listed in Ta¬ 
ble We do not include the fits for ring [d] because its 
centroid falls too close to the edge of the FOV of XMM- 
MOS at 0 ~ 14'.5 (the XMM-PN FOV extends slightly 
further but is too noisy by itself to reliably extend the 
profile outward). The inferred relative dust distances 
agree with the Chandra values to better than 5%, indi¬ 
cating that the deconvolution is robust, and consistent 
with an increase in median ring radius from Chandra 
ObsID 15801 to XMM ObsID 0729560501 of - 16%. 


As an additional consistency check, we constructed the 
radial intensity profile predicted by the deconvolution of 
the Chandra ObsID 15801 radial profile for the increased 
time delay of XMM ObsID 0729560501 (including a de¬ 
crease in surface brightness by a factor of 1.8 due to the 
^ 16% larger scattering angle, assuming a = 4). The 
predicted profiles for all four rings are plotted as dashed 
lines in Fig. 11 The profiles of rings [a]-[c] agree well 
with the XMM deconvolution, given (a) the different cov¬ 
erage of the FOV by the Chandra and VMM instruments, 
which should result in moderate differences between the 
two data sets even if they had been taken at the same 
time delay and (b) the different angular resolutions and 
the responses of the telescopes. 


3.3. CO Spectral and Image Analysis 


The deconvolution into dust sheets from §3.2| shows 
that the scattering column density must be concentrated 
in dense, well localized clouds. Gas with sufficiently high 
column density (of order 10^^ cm“^) and volume num¬ 
ber density (of order 100 cm“^) to produce the observed 
brightness and arcminute-scale variations in the echo can 
be expected to be largely in the molecular phase (e.g.lLee 
et alj2014|). 


In order to search for corresponding concentrations of 
dust and molecular gas, we compared the peaks in the 
inferred dust column density to maps of CO emission 
from the Mopra Southern Galactic Plane Survey that 
cover the entire Chandra and XMM-Newton FOV . 

While CO emission does not directly probe the dust 
distribution along the line of sight, the largest concen¬ 
trations of dust will reside in dense molecular gas, and 
CO is the most easily accessible tracer of molecular gas 
that includes three-dimensional information about the 
position of the gas, as well as a measure of the mass of 
the cold gas through the CO intensity. 

The bottom panel of Fig. shows the CO spectrum 
integrated over the entire grid of rings identified in the 
XMM-Newton image in Fig.|^ smoothed with a 2kms“^ 
FWHM Gaussian for noise removal. The spectrum shows 
five well identified, localized components of molecular 
gas in the direction of Circinus X-1. A simple multi- 
Gaussian fit (shown in color in the plot) locates the 
components at velocities ulos.i = —81kms“^, ulos ,2 = 
—73kms“^, vlos ,3 = —55kms“^, ulos ,4 = —32kms“^, 
and ulos ,5 = 9kms“^, where component 3co requires 
two Gaussian components to fit at —50kms“^ and 
—55kms“^. 
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Fig. 13.— Top panel: Radial velocity as a function of distance 
from the Sun for the sight-line in the direction of Circinus X-1 
at Galactic coordinates I = 322.12°. b == 0°. usin g the rotation 
curves from|McClure-Griffiths & Dickey|(|2007|l and Brand &; Blitz| 
( flO^ for the inner and outer Ualaxy, respectively! Uverplotted 
are the locations and radial velocities for the four dust clouds [a]- 
[d] (from top to bottom) responsible for the dust scattering rings, 
assuming the best fit distance of 9.4 kpc. Middle panel: 21cm HI 
spectrum integrated over rings [a]-[c] defined in Figl^ Component 
Ifji comprises CO components Ico a^nd 2qq listed below. Bottom 
panel: ^^CO^ectrum integrated over the same region of rings [a]- 
[c] from Fig. ^ smoothed with a 2kms“'^ FWHM Gaussian. Five 
discernible vSocity components are labeled lco-5co- 


While other components of cold gas are likely present 
below the sensitivity limit of the survey, and while some 
of the components may be superpositions of several dis¬ 
tinct clouds it is clear from the spectrum that these are 
the primary concentrations of cold gas along the line- 
of-sight and therefore the most probable locations of the 
dust clouds responsible for the observed X-ray light echo. 
We adopt the Galactic rotation curve by [McClur^ 


Griffiths & Dickey| (|2007|) for the fourth Galactic quad¬ 
rant within the boiar c ircle. We use the rotation curve by 
Brand & Blitz| (|1993|) for the outer Galaxy, matched at 
ut.sr = U following the procedure used in [Burton et al.| 
(2013), for the Galactic position of Circinus X-1. The 
resulting radial velocity curve in the direction of Circi¬ 
nus X-1 as a function of distance from the Sun is plot¬ 
ted in the top panel of Fig. Because the direction 
to Circinus X-1 crosses the smar circle, the distance- 
velocity curve is double-valued, and we cannot simply 
associate a particular velocity with a given distance. 
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Fig. 14. — Differential (black) and total (gray) Mopra spec¬ 

tra of three ring sections (see Fig. ^ for a definition of the ring 
sections), smoothed with a 3kms“~FWHM Gaussian. Differen¬ 
tial spectra show comparisons between neighboring ring sections 
that have large vs. small dust scattering intensity to isolate the 
GO velocity component responsible for the X-ray peak. Bottom: 
ring section [all] minus [al2], middle: ring section [b6] minus [b7]; 
top: sum of ring sections [c5:8] minus sum of [cl:4]. 


The rotation curve turns over at a minimum veloc¬ 
ity of Urad > —75.3kms“^. Any velocity component 
at velocities beyond the minimum must be displaced 
from the Galactic rotation curve and lie near the tan¬ 
gent point. Given the velocity dispersion of molecular 
clouds of aco ^ 5kms“^ relative to the LSR (see Q, 
velocities of Urad ^ —80kms“^, as measured for compo¬ 
nent Ico are consistent with the rotation curve and place 
the material at a distance of roughly Z? ^ 6 — 8 kpc. 

Given a rotation curve, one could attempt a global dis¬ 
tance estimate by simply comparing ring radii and GO 
velocities. However, a considerably more robust and ac¬ 
curate distance estimate can be derived by also consid¬ 
ering the angular distribution of dust inferred from the 
brightness of each of the rings as a function of azimuthal 
angle. 


3.3.1. Differential CO spectra 

In order to identify which cloud is responsible for which 
ring, we generated adaptively smoothed images of each 







































12 



15:22:00 21:00 20:00 19:00 

RA (J2000) 


Fig. 15.— Adaptively smoothed Mopra image of compo¬ 

nent 2(30 the velocity range — 76kms~^ to — 72kms“^. The 
image shows a clear lane of dense molecular gas across the im¬ 
age. Overlaid for comparison is the angular grid used to denote 
rings in Fig. Connected blue dots are the local surface bright¬ 
ness peaks in ring [all] from Chandra ObsID 15801, XMM ObsID 
0729560501, and XMM ObsID 0729560601, showing a clear spa¬ 
tial coincidence with the CO lane. Also overlaid are the centroids 
(dots) and FWHMs (arcs) of the peaks in the angular intensity dis¬ 
tributions of rings [a]-[d] from inside out in blue, green, yellow, and 
red, respectively, determined from spectral fits to Chandra ObsID 
15801 in 10-degree segments for each ring. Note that the arc for 
ring [a] denotes the secondary peak. The peak of ring [c] is close 
to the chip edge, indicated by the arrow denoting the lower limit 
to the ring extent. 


velocity component within velocity ranges centered on 
the CO velocity peaks of the clouds. The im ages are 
shown in Figs. HU and discussed in §3.3. 2| We also 
generated CO diHerence spectra to identify which veloc¬ 
ity components may be responsible for the clear bright¬ 
ness enhanceme^s seen in the different rings. They are 
plotted in Fig. which we will discuss below. 

The un-absorbed X-ray intensity shown in Figs. 1^ and 
is directly proportional to the scattering dust column 
aensity. Variations in l 2 - 5 keV as a function of azimuthal 
angle therefore probe the spatial distributions of the dust 
clouds responsible for the different rings. We can use 
a differential spectrum of two ring sections that 

have significantly different X-ray brightness to identify 
possible velocity components that show a marked 
difference in column density consistent with the dust en¬ 
hancement. 

The two most obvious regions of strong local variation 
in / 2 - 5 keV are seen at ring sections [all] and [b6], both 
of which are significantly brighter than the neighboring 
sections [al2] and [b7], respectively. 

In the bottom panel of Fig. we show the difference 
spectrum of rings sections [all] minus [al2] (i.e., 
the spectrum of [al2] subtracted from the spectrum of 
all]) as a thick black line. We also show the spectrum of 
all] only (i.e., without subtracting [al2]) as a thin gray 
line. Because of the strong enhancement in dust detected 
in [all], the corresponding CO velocity component must 
be significantly brighter in [all] as well, and therefore 
appear as a positive excess in the difference spectrum. It 


TABLE 2 

Gaussian fits to CO velocity components for rings [a]-[d] 


CO cloud 

cr[kms 

'''h 2 rir,20 f.™-2la 

2^00.20.3 


Ring [a]: 

Cloud 2(^0 

-73.56 ±0.12 1.46 ±0.12 

5.09 ±0.97 

Ring [b]: 

Cloud Ico 

-80.85 ±0.19 2.01 ±0.20 

3.88 ±0.97 

Ring [c]: 

Cloud Ico 
Cloud 2co 
Cloud 3co 
Cloud 4co 

-81.18 ±0.15 1.72 ±0.16 

-72.44 ±0.09 2.35 ±0.09 

-55.39 ±0.08 1.82 ±0.08 

-31.16 ± 0.02 2.16 ±0.02 

2.85 ± 0.35 

7.71 ± 0.40 

5.76 ± 0.34 

36.02 ± 0.38 

Ring [d]: 

Cloud Ico 
Cloud 2co 
Cloud 3co 
Cloud 4co 

-81.29 ±0.16 1.93 ±0.17 

-72.14 ±0.11 2.69 ±0.12 

-54.80 ±0.09 2.26 ±0.10 

-30.85 ±0.02 2.03 ±0.02 

3.43 ± 0.76 

7.59 ± 0.90 

7.10 ±0.78 

34.22 ± 0.73 


Note. — spectra are integrated over rings of radii 

4.0' < Ra < 6.0', 6.0' < Rb < 8.5', 8.5 < Rc < 12.0', and 
11.0' < Rd < 14.0', respectively, weighted by the 3-5 keV Chandra 
intensity. Only velocity com ponents identified in the differential CO 
spectra plotted in Fig. |14| as possible locations of dust responsible 
for the rings are listed 

hydrogen column density of the velocity component in units of 
10® cm^^/xQo 20.3i where xqq 20.3 i® the CO to H2 conversion 
factor in units of the fiducial value xco = lO^'^'^ cm“^/(K kms“t) 
l|Bolatto et al.|2013[l. 


is clear from the figure that only one velocity component 
is significantly enhanced in [all] relative to [al2], the one 
at i^rad = — 74kms“^. We can therefore unequivocally 
identify velocity component 2co at — 74kms“^ as the 
one responsible for ring [a]. 

The middle panel of Fig. 14 shows the difference 


spectrum of [b6] minus [b7], where the light echo is signif¬ 
icantly brighter in [b6] than [b7]. In this case, two 
velocity components are significantly enhanced in [b6] 
relatively to [b7], component Iqo at Urad —80kms“^ 

and component 4co at Urad ~ — 32kms“^. While we 
cannot uniquely identify which velocity component is re¬ 
sponsible for ring [b] from the spectra alone, we can rule 
out components 2co and 3co- 
Finally, ring [c] does not show a single clearly X-ray 
enhanced section relative to neighboring ones. Instead, 
the Southern sections [c5:c8] are significantly enhanced 
relative to the Northern sections [cl:c4]. We therefore 
plot the difference CO spectrum of sections [c5] through 
[c8] minus sections [cl] through [c4] (top panel of Fig. 14). 
In this case, all four velocity components are enhanced in 
[c5:c8] relative to [cl:c4], and we cannot identify even a 
sub-set that is responsible for ring [c] from the difference 
spectrum alone. 


3.3.2. CO images 

In order to compare the ring emission with the spa¬ 
tial distribution of cold gas and dust in the different 
velocity components identified in Fig. |I3[ we extracted 
images around the peak of each velocity component, 
adaptively smoothed to remove noise in velocity chan¬ 
nels of low-intensity while maintaining the full Mopra 
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Fig. 16.— Adaptively smoothed Mopra image of compo¬ 

nent Ico the velocit y ra nge —85kms~^ to —75kms“^, same 
nomenclature as in Fig. |l5| The CO emission in this band shows 
a clear peak at locations [b5] and [b6], matching the intensity dis¬ 
tribution of X-ray ring [b] (green ring segment). 


Fig. 17.— Adaptively smoothed Mopra image of compo¬ 

nent 2(30,b tti® veloci ty r ange —72kms“^ to —67kms“^, same 
nomenclature as in Fig. |15| The CO intensity peak matches the 
centroid of X-ray ring [cj (yellow arc). 


angular resolution in bright velocity channels. We em¬ 
ployed a Gaussian spatial smoothing kernel with width 
a = 1.5' X [1.0 - 0.9exp[-(z; - 'i;peak)/(2(T^o)]> where 
'f^peak and (Tco are the peak velocity and the dispersion 
of G aus sian fits to the summed CO spectra shown in 
Fig. 14 This prescription was chosen heuristically to 
produce sharp yet low-noise images of the different CO 
clouds. 

The CO intensity maps can be compared with the loca¬ 
tions of excess X-ray dust scattering to identify potential 
CO clouds responsible for the different rings, support¬ 
ing the spe ctral i dentifica tion of possible velocity com¬ 
ponents in §3.3.1 and Fig. 14 To quantify the deviation 
of the X-ray rings from axi-symmetry and for compari¬ 
son with the CO data, we constructed azimuthal inten¬ 
sity profiles of the rings from Chandra ObsID 15801 in 
ten degree bins, following the same spectral fitting pro¬ 
cedure used to construct Figs.|^and[^ Ring radii in the 
ranges [a]: 4’-6’, [b]: 6’-8’, [c]: 9’-10^5, [d]: IF.5-12’.75 
were chosen to best isolate each ring. The centroids and 
FWHM of the intensity peaks determined from Gaussian 
fits for each ring are plotted as colored circles and arcs 
in Figs. [T5]-|T^ respectively. 

It is important to note that a one-to-one match in the 
intensity distribution should not be expected on large 
scales for every ring for three reasons: (a) multiple dis¬ 
tinct clouds at different distances may fall into the same 
velocity channel because of velocity deviations of the 
clouds from the local standard of rest; (b) clouds at dif¬ 
ferent distances may fall into the same velocity chan¬ 
nel because of the double-valued nature of the distance- 
velocity curve; (c) individual rings may contain scatter¬ 
ing contributions from multiple clouds at different veloc¬ 
ities but similar distances. 

However, clear local maxima in scattering intensity 
may be expected to correspond to local maxima in CO 
emission, and we will base possible cloud-ring identifica¬ 
tion on local correspondence. Indeed, detailed matches 


exist for rings [a] an d [b], as already expected from the 
spectra discussed in [3.3.1 as well as for ring [dj: 


• Figure [T5] shows an image of CO component 2co; 
integrated over the velocity range —76kms“^ < 
■Crad < — 72kms“^, bracketing the velocity com¬ 
ponent at Urad —74kms“^ identified in the 

differential spectra as giving rise to ring [a]. 
The clear association of ring [a] with component 
2co at —74kms“^ in the CO image is strik¬ 
ing. The prominent spectral feature at u,ad 
—74kms“^ corresponds to a well-defined lane that 
runs through the position of Circinus X-1 and 
crosses ring section [all]. Overlaid as connected 
blue circles are the positions of the X-ray sur¬ 
face brightness peaks of ring [a] from Chandra 
ObsID 15801, and XMM ObsID 0729560501 and 
0729560601 (from inside out, given the increasing 
ring radius at longer time delays for later observa¬ 
tions!^ The intensity peak lies exactly on top of 
the CaJ peak and traces the CO lane as the rings 
sweep out larger radii. 


The obvious spatial coincidence and the fact that 
this is the only velocity component standing ou t in 


the differential spectrum of ring [all] in Fig. 14 


ambiguously determines that ring [a] must be pro¬ 
duced by the dust associated with the CO cloud 
2co at Urad ~ -74kms“^ 


In addition to the narrow lane in the North-Eastern 
quadrant of the image that we identify with the 
cloud responsible for ring [a], there is an additional 
concentration of CO emission in this channel in the 
Southern half of the image, which we identify as the 


Chandra ObsID 16578 is not shown because of the location of 
the peak at the chip gap and the increased noise at the chip bound¬ 
ary due to the shorter exposure time; XMM ObsID 0829560701 is 
not shown due to the overwhelming noise from the background 
flare that makes image analysis impossible. 
































cloud likely responsible for at least part of ring [c], 


and which we discuss further in Fig. 17 


This velocity channel straddles the tangent point at 
minimum velocity — 75.3kms“^ and may contain 
clouds of the distance range from 5kpc to 8 kpc (ac¬ 
counting for random motions). It is therefore plau¬ 
sible that multiple distinct clouds may contribute 
to this image and it is reasonable to associate fea¬ 
tures in this image with both rings [a] and [c]. 


Figure[l^ shows the image of cloud Ico in the 

velocity range —85kms“^ < Urad < —75kms“^, 
roughly centered on the peak at —81kms“^. The 
peak of the emission falls into sectors [b5:b6] and 
[c5:c6], while there is consistent excess CO emission 
in the Eastern part of ring [b]. This spatial coinci¬ 
dence with the observed excess scattering emission 
of ring [b] strongly suggests that cloud Ico is re¬ 
sponsible for the bulk of the X-ray scattering for 
ring [b], consistent with the difference spectrum in 
the middle panel of Fig. 14 Because of the strong 
excess foreground absorption in these ring sections 
(see Fig. 1^, a direct comparison with the intensity 
peaks in the X-ray images is more difficult than 
in the case of ring [a] (because the 3-5 keV chan¬ 
nel images are significantly more noisy than the 1-2 
keV and 2-3 keV channels, which are more strongly 
affected by absorption). However the clear X-ray 
intensity peak in sections [b5:b6] denoted by the 
green arc determined from the spectra correlates 
well with the peak in CO emission of cloud Ico- 


Figure shows a low-band image of component 
2co in the velocity range — 72kms“^ < Wrad < 
—67kms“^. The CO emission peaks in sections 
[c5:c7] and corresponds to the centroid in the X-ray 
intensity of ring [c], which is overplotted as a yel¬ 
low arc. The good spatial correspondence suggests 
that ring [c] is produced by a CO cloud at simi¬ 
lar velocity but smaller relative distance than the 


CO cloud responsible for ring [a] shown in Fig. 15 
Because of the substantial spectral overlap of both 
clouds, this image may not show the entire extent 
of the Southern part of component 2co- 


Figure [l8| shows an image of cloud 3co in the ve¬ 
locity range —60kms“^ < Wrad < —48kms“^. The 
CO emission peaks in section [c8], coincident with 
the peak in the X-rw intensity of ring [d] at the 
same location (Fig. [^. 

Given the presence of several CO components in 
the Southern half of the image, identification of the 
intensity peaks of overlapping rings [c] and [d] is not 
as straight forward. However, it is clear that not all 
of rings [c] and [d] can be explained by component 
3co alone. An additional contribution from mate¬ 
rial in other velocity channels is required to explain 
the excess emission in section [c5:c6]. Our discus¬ 
sion of Fig. suggests that the Southern compo¬ 
nent of Ico is the second component responsible 
for rings [c] and/or [d]. 



15:22:00 21:00 20:00 19:00 
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Fig. 18.— Adaptively smoothed Mopra image of compo¬ 

nent 3co iri the veloci ty r ange —60kms“^ to —48kms“^, same 
nomenclature as in Fig. |15| The CO intensity peak in the Southern 
half matches the location of the intensity peak of ring [d] (limited 
to the Chandra FOV; red arc). 

which is the dominant CO component in the Mopra 
FOV. A very bright lane of molecular gas is crossing 
the Southern half of the FOV through sections [c5], 
[b5:b7]. Comparison with the X-ray images shows 
a clear spatial overlap of this CO cloud with the 
bluest sections of the ring emission that we identi¬ 
fied as foreground absorption above. A contour of 
this cloud is overlaid on Figs. and and traces 
the region of highest column density very closely. 

Because component Aqo causes absorption in rings 

[b] , [c], and [d] without a correspondingly large 
scattering intensity peak (which would have to be 
about an order of magnitude brighter than the 
other rings, given the very large column density of 
component 4co)) the main component of cloud 4co 
must be in the foreground relative to the clouds gen¬ 
erating the scattering emission of X-ray rings [b], 

[c] , and [d]. The brightest parts of this cloud are 
almost certainly optically thick in ^^CO. The esti¬ 
mate of the cloud column density for this feature 
based on the ^^CO intensity in Table is therefore 
likely an underestimate. In fact, the differential 
hydrogen column inferred from the X-ray spectra, 
roughly TVh ~ 10^^ cm“^, is about a factor of two 
to three larger than that inferred from ^^CO. 

From the image comparisons and the differential CO 
spectra, a clear picture regarding the possible ring-cloud 
association emerges: Ring [a] must be generated by cloud 
2co at ~ —74kms“^. Ring [b] is generated by cloud Ico 
at ~ —81 kms“^. At least part of the emission of ring [c] 
or [d] is likely generated by cloud 3co; but an additional 
component is required to explain the emission in [c5:c6]. 

4. DISCUSSION 

4.1. The distance to Circinus X-1 


Figure shows an image of cloud 4co in the ve- Based on the ring-cloud association discussed in the 

locity range —33.6kms“^ < Urad < —29.6kms“^, previous section, we can use the kinematic distance es- 
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15:22:00 21:00 20:00 19:00 
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Fig. 19.— Adaptively smoothed Mopra image of com¬ 
ponent 4(30 velocity range —33.6kms~^ to —29.6kms~^, 

showing the brightest e mis sion feature in the field of view; same 
nomenclature as in Fig. |15| The band of molecular gas across the 
lower half of the im^e matches the band of absorption seen in the 
color images in Fig. and in position and morphology and corre¬ 
sponds to the excess absorption found in the spectrally determined 
absorption column shown in Figs. [^and]^ 


timate to each cloud from the velocity and the 

relative distance x to the cloud from the light echo to 
constrain the distance to Circinus X-1. 

To derive kinematic distances to the clouds, we assume 
that the molecular clouds measured in our CO spectra 
generally follow Galactic rotation, with some random dis¬ 
persion around the local standard of rest (LSR). Lit¬ 
erature estimates for the value of the ID velo city dis¬ 
persion (Tcioud range from Ccioud ~ 3.5kms“^ (Liszt & 


Burton 


(|Stark 


1981 1 for GMCs in the inner Galaxy to 7 km s 




1984 Stark & Brand 
high latitude clouds in the So 


19891 for moderate-mass 
ar Galactic neighborhood. 


Recent studies of kinematic distances to clouds assume 
a ID dispersion of aHnnri ~ 3kms“^ around the LSR 


(e.g. Roman-Duval et al. 2009 
the following, we will conserva' 


Foster et al. 2012). In 
dvely adopt a value of 


Ccioud^Skms 

In addition to random cloud motions, an unknown con¬ 
tribution from streaming motions may be present. The 
magnitude of streaming motions relative to estimates of 
the Galactic rotation curve i s generally relatively poorly 


constrained observa tionally ^ 


2009 


Reid et al. 2014). McClure-Grifhths & Dickey (2007! 
characterized the possible contribution ot streaming mo¬ 
tions to deviations from the rotation curve in the fourth 
quadrant and found systematic deviations due to stream¬ 
ing motions with a peak-to-peak amplitude of 10kms“^ 
and a standard deviation of Ustreaming 5kms“^. While 
streaming motions are not completely random, the four 
dust clouds responsible for the rings are located at suf¬ 
ficiently different distances that their streaming motions 
relative to the Galactic rotation curve may be consid¬ 
ered as un-correlated. Thus, we will treat the uncer¬ 
tainty associated with streaming motions as statistical 
and add the dispersion CTstreaming in quadrature to the 
cloud-to-cloud velocity dispersion, giving an effective ID 
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Fig. 20.— Unnormalized posterior of the distance to Circinus X-1 
from ObsID 15801 and 16578. Shown in black: Likelihood func¬ 
tion with all three priors about ring-CO-cloud associations in place 
(i.e., we identify velocity component 2qq at ulos ~ —74kms“^ 
with ring [a], component Ico 3-t ulos ~ —80kms“^ with the peak 
of ring [b], and require ring [c] and [d] to be produced by different 
CO components). The likelihood distributions without any pri¬ 
ors, with only prior one, and with prior one and two are plotted 
as dotted, dashed, and solid gray curves, respectively. The most 
likely distance of 9.4 kpc and the 1-sigma and 3-sigma uncertainties 
are plotted as solid, dashed, and dashed-triple-dotted gray vertical 
lines, respectively. 


bulk velocity dispersion relative to Galactic rotation of 
d'cloud.eff ^ 7 km s 

We assume that the mai n component of each dust 
cloud identified in Fig. and Table is associated 
with one of the GO velocity components identified above. 
However, a single velocity component may be associated 
with more than one ring, given the double-valued nature 
of the velocity curve and the possible overlap of clouds 
in velocity space due to excursions from the LSR. For 
a given a distance D to Gircinus X-1, we can then de¬ 
termine the radial velocity of each dust cloud if it were 
at the LSR, given its relative distance a;(I?, dmed), and 
compare it to the observed radial velocities of the GO 
clouds. 

For an accurate comparison, we extracted GO spectra 
that cover the FOV of Chandra ObsIDs 15801 and 16578, 
limited to the radial ranges [a]: 4.0’-6.0’, [b]: 6.0’-8.5’, [c]: 
8.5’-12’, [d]: 11.0’-14.0’ for ObsId 15801 and increase ra¬ 
dial ranges by 4% for ObsID 16578 according to eq. (10). 
We weighted the GO emission across the FOV with tne3- 
5 keV Chandra intensity in order to capture GO emission 
from the regions of brightest dust scattering. We htted 
each identified GO component with a Gaussian and list 
the relevant GO clouds used in the hts in Table [2| 

For a given assumed distance D to Gircinus X-1, we 
calculate the expected radial velocities of the cloud com¬ 
ponents and compare them to the observed GO velocities 
to Hnd the distance that best matches the observed dis¬ 
tribution of clouds. We represent each GO cloud j and 
each dust cloud i by the best-ht Gaussians from Tables 
[^and[R respectively, with dispersions crco,j and Cdust,! 
and peak positions uco,j and Udust,i, respectively. In ad¬ 
dition, we allow for an effective bulk ID velocity dis¬ 
persion of CTcioud.eff = 7kms“^ relative to the Galactic 
rotation curve as discussed above, comprised of the ran¬ 
dom cloud-to-cloud velocity dispersion and contributions 
from streaming motions. 

While the velocity distribution of each dust cloud (cal- 
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culated from the Gaussian fit to the distribution of dust 
as a function of z) is not itself Gaussian, we conserva¬ 
tively use the peak value of the velocity and assign a 
velocity dispersion (Tdist.i that is the larger of the two 
propagated z-dispersions in the positive and negative z- 
direction from the peak of the cloud. 

We combine the measurements on the ring radii from 
the deconvolutions of ObsID 15801 and 16578, corrected 


by the expansion factor of the rings given by eq. (10). 


The goodness of fit is then determined by the likelihood 
function P(D) as a function of distance D, 




... ^ exp 


( 11 ) 


(l’dust,i(71) — t’CO.ji) 


dust.i 


: + O', 


CO,j: 




cloud 


where for each dust cloud i in the product, the CO com¬ 
ponent ji with the maximum likelihood is chosen from 
the set of allowed CO components, i.e., the CO compo¬ 
nent that most closely matches the velocity t’dist.i of dust 
cloud i\ Np is a normalization factor chosen to ensure 


that / dDP{D) = 1. PjD ) is plotted in Fig. 20 
As discussed in §3.3.2| we can restrict the CO clouds 
responsible for the rings bas ed on the vi sual a nd spectral 
identification discussed in §3.3. 1| and §3.3.2| We place 
two strong priors and one weak prior on the association 
of dust and CO clouds: 

1. The first strong prior is the association of ring [a] 
with CO component 2co 

2. The second strong prior is the association of ring 

[b] with CO component Ico 

3. The third, weak prior is the requirement that ring 

[c] and [d] cannot be produced by component 3co 
alone, i.e., that either ring [c], or ring [d] are asso¬ 
ciated with a component other than 3co- 

Functionally, these priors take the form of limits on the 
possible combinati ons of rings i and clouds ji in the max¬ 
imum taken in eq. (11). That is, the first prior eliminates 


the maximum for ring i = \ and simply replaces it with 
the exponential for cloud ji = 2, while the second prior 
replaces the maximum for ring i = 2 with the exponen¬ 
tial for j 2 = 1- The third prior eliminates combinations 
in the maximum for which j = 3 for both * = 3 and 
f = 4. 

4.1.1. Distance estimate with all three priors 
With all three priors in place, eq. 0 takes the form 


P(D) = - X 


( 12 ) 


...exp 


...exp 


(t’dust,l(-P) - VCO,2) 

2 (o'dust,! + ^^00,2 + ®’cloud,eff) 

(Wust,2(-P) - VCO.lf 

2 (o'dust,2 + ®’CO,l + ®’cloud,eff) 


tt4 4 

■^i=3 maxj.^, 


(1-d 


h,3^j4,3) 


. exp 


(^dust,i(T’) - VCO,}if 


^ dust,i 


'COJi 


+ T 


cloud,eff j 


where the product of two Kronecker deltas guar¬ 

antees that rings [c] and [d] cannot both be produced by 
CO cloud J = 3, since the term in parentheses vanishes 
when j 3 = j 4 = 3. 

The likelihood distribution with all three pri ors in 
place is shown as a solid black curve in Fig. The 
68% and 99.7% confidence intervals (correspoming to 
the one- and three-sigma uncertainties in the distance to 
Circinus X-1) are plotted as well. 

The best ht distance to Circinus X-1 with 1-sigma sta¬ 
tistical uncertaintie^iE] is 

ZlcirX-i = 9.41?;® kpc (13) 

Quantities that depend on HcirX-i listed in this paper 
(such as the relative cloud distance x) are calculated us¬ 
ing the best fit value from eq. (13) and are listed in gray 
color in Table and Fig. [I^ 

The distance listed in eq. is consistent with the 
dust-cloud associations of rings [a] and [b] with clouds 
2co and Icoi respectively, it places the cloud responsi¬ 
ble for rings [c] and [d] at velocities associated with CO 
components 2co and 3co, respectively, and it places the 
CO mass of cloud ico in the foreground of all of the 
scattering screens, consistent with the observed excess 
photo-electric absorption observed in the X-ray images 
that coincides with the location of cloud ico- The loca¬ 
tions and velocities to the four scattering screens derived 
from t he b est ht distance are plotted in the top panel 
of Fig.[^ showing that the dust clouds fall close to the 
Galactic rotation curve. 

We overplot the best-ht distance and velocity scale as 
a funct ion of obser ved angle 6 in gray on the top x-axis 
Figure also shows the scattering depth 
Imed gray. Because it requires knowledge of 
x{D), this quantity is distance dependent. 

The distance derived here is both geometric and kine¬ 
matic, given that we use kinematic distances to molec¬ 
ular clouds to anchor the absolute distance scale, and 
more accurate relative geometric distances to the scat¬ 
tering screens. As such, one might refer to this method 
as either pseudo-kinematic or pseudo-geometric. Given 
that the predominant error in our method stems from 
the kinematic aspect of the distance determination, we 
will simply refer to it as a kinematic distance. 

For completeness, we will also list the constraints on 
the distance with only a sub-set of priors 1-3 in place: 

4.1.2. Distance estimate with the two strong priors only 
With both prior 1 and 2 in place, a sec onda ry peak 


in Fig. [12 
A^ndcr/ 


appears at 7.8kpc (solid gray curve in Fig. 20), with a 


Uncertainties in eq. |l3| include the effect of random cloud-to- 
cloud motions and streaming motions as discussed in |4.1| Consid¬ 
ering only random cloud-to-cloud motions would reduce the uncer¬ 
tainties to llcirX —1 = kpc, in which case an allowance for 

a systematic distance error due to streaming motions of order ±0.5 
kpc (corresponding to a velocity error of 5km/s) should be made. 
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three-sigma lower limit of Hcir > 6.3 kpc. This peak 
corresponds to the solutions for which both ring [c] and 
[d] are produced by cloud 3co- The formal one-sigma 
uncertainty of the distance to Circinus X-1 from only 
the two strong priors 1 and 2 is HcirX-i = 9.4l5'7kpc. 


4.1.3. Distance estimate with only prior 1 in place 

With only the first prior in place, the likelihood dis¬ 
tribution has two roughly equally strong peaks, one at 
7.4 kpc and one at 9.4 kpc (dashed gray curve in Fig. 20), 
and the possible distance range is restricted from 4.5 to 
11 kpc. The second strong peak at lower distance corre¬ 


sponds to a solution for which both rings [a] and 
associated with component Ico and rings [c] and 
associated with component 3co- 


b] are 
d] are 


4.1.4. Lack of distance constraints without any priors 

Without placing any priors on dust-CO cloud associa¬ 
tions, the source distance is not formally c ons trained, as 
can be seen from the dotted curve in Fig. which has 
multiple peaks at distances from 3 to 15 kpc. 


4.2. Implications of a Distance o/9.4 kpc for the 
Properties of Circinus X-1 

Our three-sigma lower limit of H > 6.1 kpc on the dis¬ 
tance to Circinus X-1 is in consistent with p revious claims 
of distances around 4 kpc (laria et al. 2005k . On the other 


hand, the 1-sigma distance range listed in eq. (13) is con¬ 
sistent with the distance range of 7.8 — 10 .5 kpcbased on 
the observat ions of type I X-ray bursts by|Jonker fc Nele-] 


mans 


(2004), which further strengthens the case for the 


kinernatic distance presented here. 

A distance of 9.4 kpc is only ma rginally larger than 
the fiducial value of 8 kpc adopted in Heinz et al. (2013). 
The age of the supernova remnant and thus the age of the 
X-ray binary increases slightly to ta w ^ 3,000y r s, wel l 
within the uncertainties estimated in Heinz et al. (2013), 
while the supernova energy would increase by about 4U7o. 

At a distance of 9.4 kpc, the X-ray binary itself 
has exceeded the Eddington luminosity for a 1.4 so¬ 
lar mass neutron star frequently during the time it has 
been monitored by both the All Sky Monitor on the 
Rossi X-ray Timing Explorer and MAXI, most recently 
during the 2013 flare that gave rise to the rings re¬ 
ported here. For example, during the peak of its bright, 
persistent outburst between about 1995 and 1997, the 
mean 2-10 keV luminosity of the source was approxi¬ 
mately Ti995.2-i0keV ~ 3 X lO^^ergss”^ (i4/9.4kpc)“^. 
We estimate the peak flux measured by MAXI dur¬ 
ing the 2013 flare to be about 2.2 Crabs, correspond¬ 
ing to a peak luminosity of Apeak,2-iokeV ~ 5.6 x 

10^®ergss“^ (£)/9 4kpc)“^. 

One of the most intriguing properties of Circin us X-1 is 
the presenc e of a powerf ul jet seen both in r adio (Stewart 
et ai.|1 99'3t and X-rays (|Heinz et al.||2007|). |Fenaer et al. 
(2004a|) reported the detection of a superluminal jet ejec¬ 
tion, launched during an X-ray flare. Taken at face value, 


the observed proper motion of /i ^ 400 mas day ^ at a 


distance of 9.4 kpc implies an apparent signal speed and 
jet Lorentz factor of F ^ /3app = Uapp/c ^ 22 (11/9.4 kpc) 
and a jet inclination of 9jet ^ 3° (9.4kpc/Zl). This would 
make the jet in Circinus X-1 even more extreme than the 


parameters presented in Fender et al. (2004a), who as¬ 
sumed a fiducial distance of 6.5 kpc. 


4.3. The Galactic Positions of Circinus X-1 and 
Scattering Screens [a]-[d] 

At a distance of around 9.4 kpc, Circinus X-1 is likely 
located inside the far-side of the Scutum-Centaurus arm. 
A location inside a (star forming) spiral arm is consistent 
with the young age of the binary. 

The CO cloud nearest to the Sun, component 4co 
at ~ —33kms“^, is likely located in the near-side of 
the Scutum-Centaurus arm, explaining the large column 
density responsible for the foreground absorption. 

CO clouds Ico through 3co are at intermediate dis¬ 
tances; association of these clouds with either the Norma 
arm (or the Scutum-Centaurus arm for cloud 3co) are 
plausible but difficult to substantiate. 

Section [TT] discussed a possible error due to streaming 
motions of CO clouds relative to the LSR as part of the 
large scale Galactic structure. For streaming motions to 
strongly affect the measured radial velocity to a partic¬ 
ular cloud, the cloud would have to be located near a 
major spiral arm. 

The dominant spiral arm towards Circinus X-1 is the 
near side of the Scutum-Centaurus arm, which crosses 
the line-of-sight at roughly 3 — 4 kpc, which we tenta¬ 
tively identified with cloud 4co- CO clouds Ico and 
2coi which dominate the likelihood distribution plotted 
Fig. are likely located far away from the center of 
that spiral arm. Thus, st reaming motions are e xpected 
to be moderate (see, e.g., Jones & Dickey 2012 and the 
artist’s impression of Galactic structure by Benjamin & 
Hurlj^) , which will limit the deviation of the observed ra¬ 
dial vSocity of cloud from Galactic rotation. Our error 
treatment is therefore conservative in including the full 
velocity dispersion of stream i ng mo tions as estimated by 
McClure-Grifhths & Dickey (2007) for the fourth quad¬ 


rant. 


Given the Galactic rotation curve plotted in Fig. fT^ 
the radial velocity corresponding to a distance of 9.4 k^ 
at the Gal actic position of C ircinus X-1 is ulsr = 


—54 kms 


-1 


Jonker et al. (2007) found a systemic radial 
velocity of Urad,Cir ~ —26±3 km s“^, which would imply a 
ID kick velocity of order Ukick.LOS ~ 25 — 35 km s“^ 


con¬ 


sistent with the values infer red for high-mass X-ray bina¬ 
ries in the SMC ( Coe|200'5 1 and theoretical expectations 
for post-supern ova center-of-mass velocities of high-mass 
X-ray binaries (Brandt & Podsiadlowski 1995). 

The potential association of Circinus X-1 with the CO 
vel ocity component 5n n at ulos ~ 9kms“^ suggested 


by Hawkes et al. 
very tar end of t 


(2014) would place Circinus X-1 at the 
re allowed distance range and require a 


very large relative velocity of that cloud with respect to 
the LSR at more than 30kms“^, both of which would 
correspond to three-sigma deviations. We therefore con¬ 
clude that Circinus X-1 and cloud 5co are unassociated. 


4.4. Constraints on the Dust Scattering Cross Section 
The most robust constraints on the scattering cross 


section, and thus the properties of the dust responsible 
for the echo, would be derived by measuring the angular 


http://solarsystem.nasa.gov/multimedia/gallery/Milky_Way_Annotated.jpi 
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TABLE 3 

Hydrogen column density for rings [a]-[d] from and 
21cm spectra for Chandra ObsID 15801 


Ring: 

AfH2 [102“cm-2] 



[a] 

S-Vi'.e 

17 1 + 17.6 

ry'y 0-1-20.3 
^'•o_io.7 

[b] 

o n+3.9 
'^•^-2.0 

13.411,7 

21.2117 

[c] 

7 7+7.7 

'•'-3.9 

oo ^+32.9 

47.9l36;| 

[d] 

7-V3.6 

20 

34.81250 


dependence of the brightness of a single ring over time. 
This is because such a measurement removes the large 
uncertainty in the total hydrogen column towards a par¬ 
ticular cloud, since it would measure the relative change 
in brightness. Unfortunately, because our observations 
lack the temporal coverage to follow the brightness of 
ea ch ring over a l arge r ange in scattering angle, as done 
by Tiengo et al. ( 2010| , we cannot probe models of the 
dust scattering cross section at the same level of detail. 

However, we can still place quantitative constraints on 
da/dn, given the scattering depth of the rings shown in 
Table and estimates of the hydrogen column density in 
both H 2 from CO and in neutral hydrogen from the 21cm 
data, we can derive constraints on the dust scattering 
cross section. The scattering angle and the scattering 
depth depend on x{9, D), and the ring-cloud associations 
depend on distance as well. 

To derive an estimate of the moleeular hydrogen col¬ 
umn density, we used the values of Nh 2 from Table [3 
adopting a value of the CO-to-H9 con version factor o' 


xco = 2 X 10^° (Bolatto et al. 20131), allowing for an 

^t 


extra factor of 2 uncertainty m xco- It should be noted 
that Xco may var y by up to an ord er of magnitude within 
indiyidual clouds (Lee et al 


IIMII ). 

al hydro 


We estimated the neutral hydrogen column density 
from the 21cm Southern Galactic Plan Survey, following 
the same extraction regions used to derive the integrated 
CO intensity: We restricted the extraction region to the 
Chandra field of view and weighted the 21cm emission by 
the Chandra 3-5 keV surface brightness. We then fitted 
Gaussian lines to the velocity components, fixing the ve¬ 
locity centroids of the components to the CO velocities 
listed in Tabled 

Because the 21cm components are much broader than 
the CO components, clouds Ico and 2co are blended 
in 21cm emission and we cannot fit them as separate 
components. To estimate the HI column density in each 
component, we divide the total column in component Ihi 
into two components with a column density ratio given 
by the ratio of the flux of the CO components Ico and 
2co- We budgeted an additional factor of 2 uncertainty 
in the neutral hydrogen column derived from the 21cm 
data. The results are listed in Tabled 

Figure shows the differential X-ray scattering cross 
section da/dVLsc per hydrogen atom as a function of scat¬ 
tering angle Sgc, evaluated at 2 keV. Despite the relatively 
large scatter, the figure shows the expected decline with 
9sc- For compariso n, the figure sh ows the scattering cross 
section derived by Draine (2003), which agrees relatively 
well with the observations in slope and normalization. 
Because we associated the entire amount of neutral hy¬ 
drogen column in the velocity range —85 to —30kms“^ 



^scatter [arcsec] 


Fig. 21.— Dust scattering cross section as a function of scat¬ 
tering angle, determined from scattering depth (Table Q and the 
column density from the measured intensity and 21cm HI 

intensity (Table plotted for Chandra ObsID 15801 (circles), 
ObsID 16578 (squares) and XMM ObsID 0729560501 (triangles) 
for clouds [a],[b],[c], and [d] in green, blue, yellow, and red, respec¬ 
tively, evaluated at energy E = 2keV. Overp lotted is the predicted 
differential cross section from |Draine| ||2003| solid black line) and 
the best fit powerlaw (dotted linej. \Ve do not include ring [d] for 
XMM ObsID 0729560501 because its centroid falls outside of the 
MOS FOV. 


to the the four clouds, the derived cross sections may be 
regarded as lower limits, which may be responsible for 
the vertical offset between model and data. 

We also plot the best-fit powerlaw fit to the cross sec¬ 
tion as a function of (dotted line): 


da 

dVt 


(5.8 zb 2.2) X 10 


-21 


2keV 


1,0007 


-3.6T0.7 


(14) 


It is clear that the large uncertainties in the total hy¬ 
drogen column dominate the scatter and uncertainty in 
this plot. Higher resolution 21cm data woud likely im¬ 
prove the accuracy of this plot significantly. Ultimately, a 
direct correlation of the scattering depth with the absorp¬ 
tion optical depth, neither of which depend on the col¬ 
umn density of hydrogen, would provide a much cleaner 
measurement of the scattering cross section. Such an 
analysis is beyond the scope of this work, however. 


4.5. Predictions for Future Light Echoes and the 
Prospects for X-ray Tomography of the ISM 
Towards Circinus X-1 

The suitable long-term lightcurve and the location in 
the Galactic plane with available high-resolution GO 
data from the Mopra Southern Galactic Plane Survey 
make Gircinus X-1 an ideal tool to study the properties 
of dust and its association with molecular gas. 

Given our preferred distance to Gircinus X-1, we can 
make clear predictions for observations of future light 
echoes. Since the four major rings described in this paper 
were produced by the molecular clouds associated with 
the GO velocity components IcO) 2cOj and 3cO) future 
X-ray observations of Gircinus X-1 light echoes will show 
similar sets of rings, with radii for a given time delay 
determined from eq. ([^ and Table 

The light echo promiced by the^O cloud with the 
largest column density (cloud 4co) was already outside 
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of the Chandra and XMMYONs at the time of the obser¬ 
vations. In order to confirm the dust-to-CO-cloud associ¬ 
ation presented in this paper, future X-ray observations 
of light echoes should include Chandra or Swift observa¬ 
tions sufficiently close in time to the X-ray flare to ob¬ 
serve the very bright light echo expected from cloud 4cO) 
which we will refer to as the putative ring [e]. Ring [e] 
will be brightest in the Southern half of the image. Be¬ 
cause Circinus X-1 exhibits sporadic flares several times 
per year, it should be possible to observe this echo with 
modest exposure times and/or for moderate X-ray flu- 
ence of the flare, given the expected large intensity of 
the echo due to the large column density of cloud 4co- 

A detection of ring [e] would improve the fidelity of 
the distance estimate primarily because it would test 
the ring-cloud-associations assumed in our three priors. 
However, the statistical improvement on the distance es¬ 
timate would be marginal because the low relative dis¬ 
tance ~ 0.25 to cloud 4co implies a large change 
in with D in such a way that D 4 ^^ = X 4 ^^D (and 
thus also the velocity constraint on cloud 4co) is rather 
insensitive to D. 

Future observations at smaller delay times will also 
provide constraints on the scattering cross section at 
smaller scattering angles. Constraints on da/dD, from 
single rings at different scattering angles are inherently 
less error-prone, as the error in the normalization of 
A^h affects all measurement equally, probing the relative 
change of scattering cross section with O^c- 

Better temporal coverage of the light echo from earlier 
times will allow detailed maps of each cloud both from 
the spatial distribution of the scattering signal and from 
the foreground absorption caused by each of the dust 
clouds, which can be compared to the scattering depth 
and the hydrogen column density from CO and HI mea¬ 
surements. 

Should a high-accuracy distance determination to 
Circinus X-1 become possible from other methods, (such 
as maser parallax to one of the CO clouds, in partic¬ 
ular cloud 4cOi which has the highest column density, 
or VLBI parallax to the source itself), the distance to 
all dust components will be determined to within 5% 
accuracy or better (set by the relative angular width 
ctheta/^^med of the dust cloud from Table [^. This would 
allow for extremely accurate tomography m the dust and 
gas distribution towards Circinus X-1 from the existing 
observations presented here and any future light echo ob¬ 
servations, providing a powerful tool to study Galactic 
structure in the fourth quadrant. 

5. CONCLUSIONS 

The giant flare and the associated X-ray light echo pre¬ 
sented in this paper are the brightest and largest set of 
light-echo rings observed from an X-ray binary to date. 
The deviation from spherical symmetry allows us to make 
an unambiguous determination of the clouds responsible 
for part of the light echo from an X-ray binary. 

Because of the full temporal coverage of the lightcurve 
of the X-ray flare by MAXI, we are able to reconstruct 
the dust distribution towards Circinus X-1 by deconvolv¬ 
ing the radial profile with the dust scattering kernel de¬ 
rived from the X-ray lightcurve. 

The kinematic distance of about 9.4 kpc we derived 
from the association of ring [a] with cloud 2co and ring 


[b1 with cloud loo i s cons istent with the estimate by 


Jonker & Nelemans (2004) based on th e luminosity of 


the t ype 1 X-ray bursts of the source (Tennant et al. 


1986). It places a large cloud of molecular gas (cloud 


Tco and the associated large absorption column in the 
Southern quadrant of the FOV) in the foreground of the 
dust screens detected in the light echo. 

At a distance of about 9.4 kpc, the source has ex¬ 
ceeded the Eddington limit for a 1.4 Mq neutron star fre¬ 
quently during the time it has been monitored by the All 
Sky Monitors on the Rossi X-ray Timing Explorer and 
MAXI, both during recent outbursts and, more impor¬ 
tantly, in average luminosity during its persistent high- 
state in the 1990s. 

T aken at face value, th e observed pattern speeds in the 
jet ( Tender et al.||2004b ) indicate a very fast jet speed of 
F > 22 and a correspondingly small viewing angle of 

0jer^ 2.6°. 

Future measurements of light echoes from Circinus X-1 
will allow accurate measurements of the dust scattering 
cross section as a function of scattering angle and an 
estimation of the dust-to-gas ratio in the different clouds 
identified in the field of view towards Circinus X-1. 

The technique used to derive the kinematic distance to 
Circinus X-1 holds significant promise for future obser¬ 
vations of light echoes from X-ray transients. 
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